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COMPUTER PROGRAM FOR CALCULATING PRESSURE-BROADENED 
RAMAN SPECTRA FOR MOLECULAR NITROGEN AND OXYGEN 


by Gustave C. Fralick 
Lewis Research Center 

SUMMARY 

A computer program is given for calculating the rotational Raman spectra for mo- 
lecxilar nitrogen and oxygen. Provision is made for pressure broadening. Several sam- 
ple results at various pressures and temperatures are shown. The relative heights of 
some of the lines are affected by pressure broadening, and in nitrogen the odd-j lines 
are most affected. For instance, the j = 5 anti-Stokes line increases about 15 percent 
relative to the j = 6 Stokes line for a change in pressure from 1 to 40 atmospheres, 
while the j = 6 anti -Stokes line shows no change. 


INTRODUCTION 

In the development of gas turbine engines, there has been a continuing goal to further 
increase turbine inlet temperature and thereby increase engine performance. The mea- 
surement of gas temperatures in the region of the turbine inlet and combustor exit is re- 
quired to see that the combustor generates the right temperature profile as well as tem- 
perature level. Traditionally, these temperature measurements have been made with 
thermocouples. Presently, however, no thermocouple is available which will withstand 
the high temperatures and pressures (up to 2400 K and 40 atm) predicted for future tur- 
bine and comb\istor development. 

Other methods, such as cooled immersion probes and gas analysis, have been con- 
sidered as alternative solutions for making high gas temperature measurements, but as 
yet none have been used on a routine basis. Another method of cimrent interest is a 
remote-sensing light-scattering technique utilizing Raman scattering. The Raman scat- 
terii^ technique is attractive in that it can be used to focus remotely on a small volume 
in a gas stream and extract information associated with local gas density, composition, 
and temperature. 


Over the last few years several groups have been conducting development work on 
this technique (refs. 1 to 5), and, in at least one case, the technique has been applied to 
the measurement of jet engine exhaust temperatures up to 955 K (ref. 6). 

Some of the problems associated with the use of Raman scattering in measuring 
temperature in jet engine exhaust gases are fluorescence, ambient light background, in- 
terference of various species of gas, and pressure broadening. 

This investigation was conducted to extend the results found in reference 2 to pro- 
duce a computer program for computing the Raman spectra which includes the effects of 
pressure broadenii^. This program can be used to predict the extent of pressure broad- 
ening associated with temperature measurements at elevated pressures. 


THEORY 

The theory of rotational Raman scattering and line broadening is discussed in ref- 
erences 1 and 2, so that only a relatively brief outline is given in this section. 

When a beam of light, from a laser, for example, is scattered from a molecular 
system, most of the light is scattered with no change in wave number (the wave number 
is the inverse of the wavelei^th). This is called Rayleigh scattering. Some of the scat- 
tering molecules, however, are induced into undergoing a change of energy level by the 
light beam, and the scattered light has a wave number different from that of the incident 
light. This is called Raman scattering, li the molecule undergoes a transition to an 
energy state lower than it originally had, the wave number of the scattered light is in- 
creased (anti-Stokes scattering); conversely, if the molecule undergoes a transition to a 
higher energy state, the wave number of the scattered light is decreased (Stokes scatter- 
ing). Because the molecular energy levels are quantized (discrete), the wave number 
shifts are also discrete, and the spectrum of the scattered light consists of the unshifted 
Rayleigh line and a series of lines on either side of it, the Raman lines. The magnitude 
of the frequency shifts and the relative intensities of the various lines depend on what 
type of molecule is involved and on the temperature. 

The scattered light is not observed as a series of perfectly sharp lines, but as lines 
which are spread out, or broadened. The broadening has a variety of causes, among 
which are the fact that neither the light source nor the detector is perfectly monochro- 
matic and the fact that the molecules are in motion rather than being at rest. Since the 
broadening depends on both the temperatiire and pressure of the gas, it must be taken in- 
to account if Raman scattering is to be used as a tool for measuring the temperature of 
high-pressure gases. 

K the Raman scattered light from a molecule has an intensity I-, where j is the 
quantum number associated with a particular line, the broadening of that line is described 
by writing 
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= I^S{v, i/j) 


( 1 ) 


That is, the intensity at wave number v due to the line located at ia and having an in- 
tensity I. is found by multiplying I. by the factor S(v,v^). This is shown in figure 1. 
The factor 8(1 ^, ia) is called the normalized shape function (ref. 2). (Symbols are defined 
in appendix A. ) 

For a series of spectral lines, the intensity at v is just the sum of the intensities 
due to each line 


1(v)=TjIM=YjI.S(v,u.) ( 2 ) 

j J j J J 

The Raman spectrum is the intensity at each value of v, so the problem of calculating 
the spectrum reduces to calculating the line intensities Ij and the normalized shape 
function S(i/, ia). The shape function S(i^, i/j) is normalized to unit area in order to sat- 
isfy the requirement of conservation of energy. The total intensity of the broadened lines 
over all wave numbers is 



8 ( 1 ^, v^)dv = 



h 


which is the total intensity of all the unbroadened lines. 

The Raman lines on the low -wave -number side of a Rayleigh line are called the 
Stokes lines, and for the rotational scattering from the diatomic nitrogen and oxygen 
molecules (N2 and Og) considered in this report the line intensities are given by refer- 
ences 2 and 7 as 


l(S) ^ c (3 + l)(j + 2) e 
J ^ 2j + 3 T ^ ^ 


j - 0, 1,2,3, . 


(3) 


The lines oh the other side of the Rayleigh line are called the anti -Stokes lines, and their 
intensities are given by 


(AS) ^ c. §. e-^<j+2)(j+3)/T 

T ^ 2j + 3 T ^ 


j = 0,1, 2, 3... 


(4) 


where Cj contains other factors, such as nuclear spin. The spin causes the intensity 
of every other Ng line to be reduced by 50 percent and causes every other O2 line to be 
absent. The wave number Vq is the wave number of the laser, and 
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( 5 ) 


^ he V 2 


is the wave number shift. The lines are located at 


V. = Vq ± Av. 


( 6 ) 


The constant 0 is called the rotational characteristic temperatime of the molecule. 

Line broadening is due to the Doppler profile in the laser , combined natural and 
pressure broadening in the molecular gas, and instrument distortion (ref. 2). The pro- 
file of the laser line is given by a Gaussian function (refs. 8 and 9) 


Sl(v, i/q) = 








(7) 


and the natural -pressure broadening by a Lorentz profile (ref. 10). 


R 


R 


Sr( V, 1/’) = — 


V 2 / 


( 8 ) 


The combined effect of the Doppler profile and natural -pressure broadening is given 
by the Voight broadenii^ function, which is the convolution of and Sr: 




(9) 


Note that the two broadening functions Sj^ and Sr are normalized to unit area. For 
this reason, Sy is also normalized to unit area: 
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The two numbers and Rj^ are called, respectively, the laser and Raman band- 
widths and give the width of a profile when the he^ht is reduced to one -half of the peak 
value, as shown in figures 2 and 3. The bandwidth R^ is constant, but the bandwidth 
% is not. It varies with temperature, pressure, and the quantum number j. The form 
of the variation, deduced from the data in references 11 and 12, is 

( 10 ) 

where Tq is a reference temperature, T is the gas temperature, and P is the gas 
pressure. The constants A, B, and Tg are fitted to experimental data meas\ured by 
using gases of known temperature and pressure. 

The behavior of line broadening due to the detecting instrument is given by the slit 
function (ref. 13) 




(jA + B)P 



V^Rg 

Finally, the combined effect of all the broadening is the convolution of the Voight and slit 
functions to yield the overall shape function 
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The shape function is also normalized to unit area, as explained previously. 

In principle, the spectrum can now be calculated by substituting I. from equation (3) 
or (4) and S{v, v^) from equation (12) into equation (2). Unfortunately, the integrals in 
equation (12) cannot be evaluated analytically, but must be evaluated numerically. 

The inner integral in equation (12), the Voight integral, can be put into another form, 
for which mmierical solutions now exist. Let 



With the change of variables 



The result is 



(13) 


(14) 


(15) 


Now let 


R 

0 = 2ff;3 = 2\^n~2 — 

Rt 


(16) 
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The Voight integral may be expressed in terms of the first Doppler broadening in- 
tegral (ref. 14) 

ia) = I i//(x, 0)| = -L. p(x,e) (17) 

^ ^ L 0 J 

and aside from a constant factor independent of j , the shape function becomes 



The j -independent part of the normalization constant is dropped because only the relative 
peak heights are of concern. 

Since reference 14 supplies a computer program for numerically evaluating i//(x, 0) 
(listed in appendixes B and C as PSIPHI), the shape function can now be found numer- 
ically as soon as realistic limits of integration are foimd to replace the infinite limits 
shown in equation (18). The details of finding these limits are given in appendix D. 


DISCUSSION OF RESULTS 


Eight sample results are included in figures 4 to 11. These are results for Ng and 
Og at 1, 20, and 40 atmospheres and 300 K, for Ng at 1 atmosphere and 1100 K, and for 
Ng at 19.4 atmospheres and 300 K. The results are normalized so that the highest peak 
has the value of 10. The laser line was taken as 5145X10"^^ meter (5145 A) (Ar"*" 
laser), and the spectra were calculated on a 30x10" -meter (30-A) interval on either 

1 A o 

side of the laser line (60X10 m (60 A) for Ng at 1 atm and 1100 K). Note that the cal- 
culated spectra do not include the Rayleigh line. The following values were used for the 
constants: 


Molecular rotational characteristic temperature 0 , K 

For nitrogen 

For oxygen 

Laser bandwidth, Rj^, cm“^ 

Slit bandwidth, Rg, cm“^ 

Pressure-broadening constant. A, cm"^/atm 

Pressure-broadening constant, B, cm“^/atm 

Reference temperature for Rj^, Tq, K 


. 2.862 
. 2.065 
. 0. 166 
. 3.31 

-0. 0025 
+0.091 
. . 303 
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The values for A and B, the pressxire -broadening constants, and Tq, the reference 
temperature, were obtained from the data given in references 11 and 12. The same 
values were used for both nitrogen and oxygen. A Raman spectrvim for Ng at 19.4 at- 
mospheres and 300 K measured in a high-pressure furnace at the Lewis Research Center 
is included in figure 11 for comparison. The agreement between it and the calculated 
spectrum shown in f^ure 11 is satisfactory, although better values for A and B have 
yet to be determined. The fact that the locations of the higher j lines are different in 
the two spectra indicates that there is also a small error in e. 

As can be seen, the effect of pressure broadening is considerable, at least on the 
odd-i lines. For instance, the j = 5 anti-Stokes line in Ng increases in relative height 
about 15 percent as the gas pressure is increased from 1 to 40 atmospheres. The j = 6 
anti-Stokes line does not increase at all. The odd-j lines are more affected by broad- 
ening because they are between two tall lines, and the even-j lines are between two 
short lines. 

Increasing the temperatmre has the effect of increasing the height of all the lines and 
shifting the highest intensity line toward higher j values. At 300 K, the highest line is 
the i = 6 Stokes line, and at 1100 K the highest line is the j = 12 Stokes line. 


CONCLUDING REMARKS 

A computer program is presented which is capable of calculating the rotational 
Raman spectrum of either N 2 or O 2 and which includes the effect of pressure broadening. 
Althoxigh a general spectrum can be computed, the best method of using these calcula- 
tions to determine the temperature of a high-pressure gas has not yet been determined. 
Nevertheless, the sample calculations (figs. 4 to 11) show that the effects of pressure 
broadening are not negligible for the pressures of interest in future combustor applica- 
tions (up to 40 atm). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 11, 1976, 

505-04. 
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APPENDIX A 


SYMBOLS 

pressure -broadening constants, cm“Vatm 
spin factor for unbroadened lines 
speed of light 
Plank constant 

intensity of broadened line at wave number v 
intensity of unbroadened line, (W/cm )/cm“ 
intensity of anti -Stokes line 
intensity of Stokes line 

intensity associated with broadening of single line 

rotational quantum number 

Boltzmann constant 

pressure, atm 

laser bandwidth, cm"^ 

natural -pres sure bandwidth 

slit bandwidth 

) normalized broadening function 
laser broadening function 
natural -pressure broadening function 
slit broadening function 
Voight function 
temperature, K 
reference temperature for Rjj 
parameter in Doppler broadening integral 
parameter in Doppler broadening integral 
molecular rotational characteristic temperature, K 
wavelength, 0. 1 nm (A) 
wave number, cm"^ 


Aia 

*^0 

wave number 

wave number 

wave number 

0) 

first Doppler 


of iihbroadened line 
shift of line 
of laser line 
broadening integral 
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APPENDIX B 


DESCRIPTION OF PROGRAM 

The computer program is written in FORTRAN IV for use on the IBM 360. 

The main program is called MAINl. In this program, all the necessary constants 
are read in, and all the various parameters needed by the subprograms are calculated, 
as are the number and location of the rotational lines for either oxygen or nitrogen. The 
program uses wave numbers (cm"^) in the calculations and determines the intensity 
Mv) = ia) at each value of v. The value of L is calculated by the subroutine 

XNITRO. ^ 

Equation (AlO) indicates that the bandwidth of S(v, v^) is roughly fRy + Rg, so com- 
puting time is re dyed by excluding any value of ia (the center of the line) which is 
fartherihan 3, + r| from V from the above sum. For convenience in plotting, the 

output is in terms of I(A.), where X is the wavelength in tenths of a nanometer (ang- 
stroms). 

The actual calculation of S( v, i/j) takes place in the function subprogram called 
SHPFN. There the integration limits are calculated for each value of u. SHPFN calls 
the function PROD, which in turn calls the function VOIGHT, which uses subroutine 
PSIPHI to produce V^x, 0). PROD also uses v the function SLTFN to produce Sg(i/V). 
The function PROD is the product of Sg and if/. SHPFN uses the integration routine 
QUADl to integrate PROD between the limits uq ± 3By^. The result is then SHPFN = 

S(i/, i/j). The slit function was made a separate subprogram so that it could easily be 
changed if desired. For instance, a Lorentz distribution instead of a Gaussian function 
could be used. 
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APPENDIX C 


PROGRAM LISTING 


c... H&IHI^HAIN PROGRAn FOR CALCULATING. ROTAtlORAL RARAN SPECTRUR OF N2 OR 02 
CORflOK/BA SIC /CENTER 
COHHON/SPEC/BATIO,IHALP, LEH, IBD 
COHHON/BAR/RR.RLfRS 
CORBON/VPAFAR/THET 
COHHOR/SLPB AH/CON 

COHHON/SFARAH/RVSQ, RSSQC, OENOH.TRBBEB 

DIRENSION RSBR (100) ,BS8VSQ(100) ,THRB (100) ,DEBR(100) , CTOPP ( 100) ,THT (100) 

DIHBNSION DEINO (100) , DNO (2) 

DIRENSION ANG(IOOO) ,TENS (-1000) ,COUNT(1000) 

INDEXsl 

1 HBITB (6,101). 

101 P6RRAT(1H0» STEP 1-ROTATIONAL RARAN SPECTRUR OF NITROGEN (ENTER 1,) OR OXTGEN (ENTER 2.)*) 

BEAD(5,201) EL 

LEH-iPIX (EL) 

2 NRITE(6,102) 

102 PORHAT (1 HO'STEP 2-NHAT IS. THE ROTATIONAL CHARACTERISTIC TERP.(DEG. K)?*) 

READ(5,201) THETA 

201 FOBHAT (P10.2) 

C1»2.70O0*THBTA 

3 WRITE (6,103) 

103 PORHAT(1HO* STEP 3-ENTEB THE LASER AND SLIT BAND9IDTHS (INV. CH.)*) 

BEAD(5,201) RL 

READ(5,201) RS 
RSSQsRS^RS 
C0M*2. 7725087/BSSC 
THB*1.6651092/RL 
U HRiTE(6,10U) 

104 FORRAT(1HO*STBP 4-HHAT IS THE WAVELENGTH OF THE. LASER LINE (ANGSTRORS) ? • ) 

READ(S,201) ALAHO 

5 WRITE (6^1 06) 

105 POFHAT(1H0»STEP 5-OVER WHAT RANGE (ANGSTRORS) ON EITHER SIDE OF THE CENTER IS THE INTENSITY TO BE CALCOUTED?') 
READ(5,201) ALRRNG 

ALHRIN=ALAHO-ALHBNG 
AL HHAX= ALA HO ALRRNG 
XNOO='1.E*08/ALAHO 
XNUNAX=1 . E408/ALHHIN 
XN0HIN>1.E+08/ALRHAX 
XNOBNG«XNOHAX • XNUO 

C... CALCULATE THE HDHBER OF ROTATIONAL LINES 
BJ=AINT (XNUHNG/CI - 1.) 

BJ=B0^1. 

NTEBRsIFIX (BJ) 
e WRITE(6, 106) 

106 FOBHAT (1H0 *STBP 6-AT HOW MANY POINTS IN THS RANGE IS THE INTENSITY TO BE C AtCOLATED ( 500. RAX.)?*) 

READ(5,201) PTNO 

NPTNO=IFIX(PTNO) 

DNU (1)sXNUBNG/PlNO 
DNa(2) »(XMUO-XNQHIN)/PTNO 
7 WHITE(6,107) 

107 POFHAT(1HO STEP 7-EHTER THE GAS TEHP. (DBG. K) AND THE GAS PRESSOR E (ATH .)• ) 

READ(S,201) TEHP 

READ (5,201) PRESS 
TPPAC*SQHT(303./TEHP) *PBESS 
RATI 0*THBTA/T EBP 

6 WRITE(6,1O0) 

100 FOBHAT (1H0*STBP 8-ENTEfi THE PRESSURE BROADENING CONSTANTS, A $ BM 


READ(5,201) k 
BEAD(S«201) B 

C... CALCOLATE VABIAELES BUICK DEPEND ON MNE NUHBER 
DO 1000 KL=l,NTEBn 
K=KL-1 
AK=FLOAT (K) 

C... CALCULATE PRESSURE BROADENED BAHAK BANDBXDTH 
RR=TPFAC*(A*AK ♦ B) 

RSBB (KL)=BR 
RB02=RR/2. 

C... CALLULATB VOIGBT BABDMIDTH 

RV*RR02 ♦ SQBT(RR02**2 ♦ BL**2) 

RVSQ=RV**2 
RSBVSQ (KL)sRVSQ 

THREEB*3./SQFT(1,/RSSQ ♦ l./SVSQ) 

THRB (KL)=THBEEB 
DENOM=RSSQ-. ♦ R VSO 
•- DENfl (Kl)-DBNOJl 
CTOFP(KL) s3. ♦SQRT(DEHON) 

C... CALCULATE VOIGHT PABAHETEfi 
THET»THE*RR 
THT (KL)=tRFT 
DELNU (KL) *C 1* (AK ♦ 1 . 5) 

CESTER=XHUO ♦ DELHU(KL) 

RSSQC*RSSQ*CENTER 
1000 CONTINUE 

C... CALCULATE SPECTRUM. IHALP = 1 IS ANTI-STOKES SIDE 
DO UOOO IHALF « 1,2 
SIGN*1. 

IF (IHALF .EQ. 2) SIGHs-1. 

DO 3000 H*1,NPTNO 

XNUU«XNUO ♦ SI6N*FL0AT(H) •ONU(IHALP) 

L*NPTNO^M+1 

IF (IHALF .EQ. 2) L=NPTNO ♦ H 
ANG(L) *1.E*08/XN0U 

sunso. 

DO 2000 Is1,NTBnfl 

J«I-1 

INOa2 

IF(2*(J/2) ,NE. J) IND»1 
C... INOal IP J IS ODD,INOa2 IF J IS' EVEN 
XJapLCAT (J) 

RRsRSBR (I) 

RVSQa RSBVSQ (I) 

THREEBaTURB (I) 

DENOnaDENH (I) 

CUTOFF=CTOFF (I) 

THETsTHT (I) 

CENTERS XNUO ♦ SIGN*DELN U (I ) 

RSSQCaRSSQ*CENTER 
TSTaABS(XNUU - CENTER) 

C... SKIP LINES WHICH ARE TOO FAB AWAT TO CONTRIBUTE 6 SKIP EVEN 02 LINES 

IF ( (TST ,GE. CUTOFF) .OR, (LEM . EQ . 2 .AND. IND .EQ. 2) ) GO TO 2000 
CALL XNITRO (XJ,PK) 

SUM a sun ♦ PK*SHFFN(XNOU)/BB 
2000 CONTINUE 

TENS(L) =SUH 
3000 CONTINUE 
4000 CONTINUE 
9 WHITE(7,110) INDEX 

IF(INDEX .NE. 1) GO TO 11 

c... iNDEXa 2 Tells plotting routine no hore data is to follow 

c... NORMALIZE 6 PEED INTO DATASET "DATA1" FOB PLOTTIBG,OSmG WRITE (7, . , • ) 
NTa4*NPTNO 

C... HT IS TOTAL NO. OP POINTS,COONTS S WAVELENGTHS, FROM BOTH SIDES OF THE SPECTRUM 
NTOTs2*NPTNO 

C... FIND LARGEST VALUE 
AMAX=TENS(1) 

DO 5000 11=2, NTOT 

IF(TBNS(II) .GT. AMAX) A BAX=TE NS (I I ) 

5000 CONTINUE 

C. .. SCALE SO THAT LARGEST VALUE IS 10 
XNCBM=i 1*AMAX 
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110 

111 

c. .. 
c. . . 

6000 

109 

202 

10 

11 

c. . . 
c« , . 

1 

c. .. 

2 

c. . . 
c... 


c. • . 


1 

2 

C. . . 


HBITE(7,110) LEM 
HBITE (7,110) NT 
PORHAT (14) 

WRITE (7,111) TEMP 
MRITE(7,111) PRESS 
PORHAT (E13.6) 

DO 6000 K=1,NTOT 
COUNT (K)=TENS (K)/XN0BH 

DATA IS PED IN SERIALLY. THE PIRST NO. IS THE HEIGHT THE SPECTROM AT POINT K,COOBT (K) ,AND THE 
SECOND IS THE LCCATtON OF THE POINT IN ANGSTROBS,ANG (K) . 

HRITE(7, 111) COONT (K) 

BRITE(7,111) ANG(K) 

WRITE (6,109) 

PORHAT ( 1HO»TO WHICH STEP DO TOO WISH TO RET0RH7-9 GETS YOU OOT») 

READ (5,202) HL 
PORMAT(II) 

GO TO (1, 2, 3, 4, 5, 6, 7, 8, 10), HL 

INDEX*2 

GO TO 9 

STOP 

END 


SOBROUTINE XHITRO (AO , PKH T) 

CALCULATES PEAK HEIGHT, PKffP, POR H ITROGEN (LEH= 1 ) 6 POR OXYGEN (LBR* 2) FOR ODD VALDES OP J 
COHHON/SPEC/RATIO,IHALP, LEM,IND 
COHHON/BASIC/CENTER 
OVERALL CONSTANT CHITTED 

C=RAT10*(AJ*U) * (AJ*2. 1 * (CENTER»*lt)/(2. *AJ * t 
IP (IHALP . EQ. 1) GO TO 1 
PKHT«C*EXP (-RATI0*AJ* (AJ + 1.) ) 

IF(LEH .EQ. 2) RETURN 
GO TO 2 

PKHT=C*EXP (-RATIO* (AJ*2 . ) • (AJ*3. ) ) 

IP (LEH . EQ. 2) RETURN 

ODD 0 LINES ARB HALF INTENSITY POR NITROGEN 

IF(IND .EQ. 1) PKHls0.5*PKHT 

RETURN 

END 


FUNCTION SHPFN (XNU) 

CALCULATES UNNORHALIZEO SHAPE FUNCTION AT WAVE NO. X NO 

COBHCN/SPARAH/BVSQ, RSSQC, DENOH,THBEEB 

COHHON/WVNO/GNU 

EXTERNAL PROD 

GNU-XNU 

CALCULATE INTEGHATION LIHITS 
CENT=(BVSQ*XND ♦ RSSQC)/DENOH 
XHIN=CENT - THREEB 
XHAX=CENT ♦ THREEB 
SHPFN=Q0AD1 (5 , 3 , X BIN, X H AX , PROD) 

RETURN 

END 


FUNCTION PROD (XNP) 

PRODUCT FUNCTION NEEDED BY QUADi 

CCHHON/FLAG/KF 

COHHON/WVNO/GNU 

XNNsGNU 

FAC1*SLTFN (XNP,XNN) 

IF (KF . EQ. 2) GC TO 1 
PROD=PAC1*VOrGHT(XNp) 

GO TO 2 
PROD=0. 

BETURN 

END 


FUNCTION VOIGHT (XN) 

CALCULATES VOIGHT FUNCTION BY USING PIRST DOPPLER BROADENING INTEGRAL 


COnflCN/BAfl/Ra,BL,flS 
COBPION/BASIC/CENTHB 
COnnON/VPABAM/THET 
C... OVERALL CONSTANT CHITTED 
X=2.*(CENTEB-XN)/RR 
CALL PSIPHI (X,THET, PS I, PHI) 
VOIGHTsPSI 
RETURN 
END 


FUNCTION SLTPN (X,XS) 

C... GAUSSIAN SLIT FUNCTION CENTERED AT XS, BANDWIDTH RS 
COHHON/SLPBA51/CCN 
COMHON/FLAG/KF 
ARG=CON* { (X-XS)**2 ) 

IF(ARG .GE. 170.) GO TO 1 
SLTPN=EXP (-ARG) 

KP-=1 
GC TO 2 

1 SLTFN=0. 

KF=2 

2 RETURN 
END 


SUBROUTINE P SI PH L (X X , TT , U , V) 

C... REF; B. CANPIGHT C T. SEHLER, TH X-2559 

C... XX s 2.*(NU SUB J - HU)/RR 

C... TT = THETA = 2.*SQRT( LN(2.) )*HR/RL 

C... U r. V ARE FIRST AND SECOND DOPPLER SROADBNtNG INTEGRALS 

C, . . 

C... H«1. 

OIHENSION EHN(S) ,EN2(10) 

DATA PI22/12. 566371/ 

DATA RTP r/1 . 77 245 39/, PC H /I .5707 963/ 

DATA PI,PI2,N/3. 1415927,6.2H31854,5/ 

DATA ERN/0., . 11709966E0, . 583004 89 E- 2 3928256 1 E-4 , - 

1 . 35821059E-7/ 

DATA EN2/0. , .24789999, .1 1 709966 ,. 3354 96 1 5B-1 58300489 E-2, - 

2 . 6 14U82 64 E-3, . 39 282561 E-4, . 1 5 2 3 1 50 2B- 5, . 35820159E-7, . 51 095996E-9/ 

XDaAbSlXX) 

THETAaADS (TT) 

YsTHETA/2. 

X»XD*Y 

C... ARE HE IN DANGER? 

XTEST=X-AINT (X) 

IF( (XTEST .LT. .01 .OR, XTEST .GT, .99) .AND. ABS<Y) ,LT. .01) GO TO 40 
IP( XD .GT. 75. *SORT(1 , ♦ 1 ./ (Y*Y) ) ) GO TO 30 
IF (Y .LT. .001) GC TO 40 
C... BEGIN SERIES SUHHATION 
Y2 = Y*Y 

C. . . SET N=0 TBFH 

D=FI* (X*X + Y2) 

SU=Y/D 
SV=X/D 
DO 10 1=2, N 
AH=FLOAT (1-1) 

XNP=X-AH 
XNH=X*An 
DP=XNP*XNP ♦ Y2 
DB=XNR*XNH ♦ Y2 

SO = SO ♦ EBN (I) *Y* (1 ,/DP ♦ 1./DH) 

SV = SV ♦ EBN (I) ♦ (XNP/DP ♦ XNM/DH) 

10 CONTINUE 

IP (Y .GT, PI) GC TC 20 
P = 2. 

IF (Y . BQ. PI) P=1. 

C... BEGIN POLE CONTRIBUTION 
XP2=X* PI2 
ARG1=Y*PI2 

IP(ARG1 .LB. -170.) GO TO 14 
EYP2 = EXP (ARG 1) 


A 


15 


o n n o o n 


GO TO 15 
ia EYP2-0'. 

15 COHTIHUB 
XY2*X*Y*2. 

SXY2>SIN (XY2) 

CXY2«COS(XY2) 

SXP2sSIH (XP2) 

CXP2*COS (XP2) 

ARG2*Y2 - X*X 

IP(ABG2 .LE. -170.) GOTO 16 
EYX*BXP(ARG2) 

GO TO 17 

16 EYX=0. ' X 

17 0*1.- BYP2* (2. *CXP2 - EYP2) 

SO * SO ♦ P*BYX*( CXY2 - BYP2»(CXP2«CXY2 ♦ SXP2*SXY2) ) /O 
SV s SV - P*EYX*( SXY2 ♦ EYP2* (SXP2 ♦CXY2. CX^2*SXY2) )'/0 
20 0*S0*Y«RTPI ' 

V*SV*THBTA*RTPI 
IP (XX .LT. 0.) V*-V 
60 TO 50 
30 paRTPI*Y 

SO a 1./( (1. ♦ -XD*XD)*P ) 

SV=S0*XD 
60 TO 20 
40 CONTINUE 

IP(XD .LT. 4. .AND. THETA .LB. .009} RRITE(6« 101) X D^TRETA 
101 PORHAT (IHO’XD * • El 3. 7, * T HET A * *E13.7,*PR0H PSIPHI*) 

C... THIS IS EXPANSION NITH H*0.5 
Y2*Y*Y 

D*PI» (X*X ♦ Y2) 

SO*Y/D 
SVa x/D 

DO 45 I»2, 10 
AH*0.5*PLOAT (1-1) 

XNPaX-AH 
XNfiaX^AH 
DP*XNP*XHP ♦ Y2 
DaaXNn«XHH ♦ Y2 

SO a SO ♦ EN2(I) *Y* (1 ./DP ♦ 1./DR) 

SV a SV ♦ EM2(I) *(XNP/DP ♦ XNa/DH) 

45 CONTINUE 
SU»S0*0.5 
SV»S?*0.5 

IP (Y .GT. PX2) GO TO 20 
P*2. 

IP (Y .EO. PI2) E=1. 

XP2aX*Pl22 

ABG3*Y*PI22 

IP(ARG3 .LE, -170.) GO TO 40 
EYP2*EXP (ARG3) 

GO TO 15 
48 EYP2*0. 

GO TO 15 
SO RETURN 

END 


FUNCTION Q0AD1 (NN,NS, XI , XOrPOPX) 

NOnERICAL INTEGRATION USING THE HETBOO OP GAUSSIAN QUADRATURE 
NN=NO. OP POINTS/DIVISION, 00*3, 4, 5, . . . , 16 

HS*NO. OP DIVISIONS INTO BHICH THE INTERVAL (XO - XI) IS DIVIDED 
XI=LOVBR LlaiT CP INTEGRATION 
XO=DPPEB LIHIT OF INTEGRATION 

POFXaUSEB SUPPLIED EXTERNAL FUNCTION SOBPROGRAR 
DIHENSION C ( 140) ,A (70) , H (70) , A A (56 ) , B6 (56) ,CC(28) 

REAL INTGRL 

EQUIVALENCE (AA (1 } ,C ( 1) ) , (B6 (1 ) ,C (57 ) ) , (CC (1) ,C (1 13) ) 

DATA AA/- 

1 7.74596669241483E-01, 5 . 55555555555556 B-01 ,-0 . B-00,- 

1 8.08888888808889E-O1, 8 . 6 1 1 363 1 1 594053E-01, 3, 47854045137454B-O1 
1 3.39981O43584056E-O1, 6 . 521 451 54 862546E-01 , 9.06 179845938664B-01, - 
1 2.36926885056189E-01, 5. 384693 10 10 5683E-01, 4.7062867O499366B-O1 
1-0. E-OO, 5.68888088888889E-O1, 9. 32469514203 152E-01-- 
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o u 


\ 1.7132ft492379l70E-01, 6 . 61209396466265E-01, 3.60761573048139E-01 
1 2.38619186083197B-01, 4 .67913934572691E-01 , 9. 49107912302759E-01,- 
1 1.2908a966168870EH)1, 7 . 4 1 S3 11 S5599394E>0 1, 2.797053914892778-01,' 
1 4.058051513773978-01, 3.818300505051198-01,-0. BrOO,- 

1 4. 179591836734698-01, 9.602898564975368-01, 1.012285362903768-01,- 
1 7.966664774136278-01, 2.223810344533748-01, 5.255324099163298-01,- 
1 3.137066458778878-01, 1.834346424956508-01, 3. 626837033783628-01 ,- 
1 9.681602395076268-01, 8.127438836157408-02, 8.360311073266368-01,- 
1 1.806401606948578-01, 6.133714327005908-01, 2.606106964029358-01,- 
1 3.242534234038098-01, 3.123470770400038-01,-0. 8-00,- 

1 3.302393550012608-01, 9.739065285171728-01, 6.667134430868608-02,- 
1 8.650633666889858-01, 1.494513491505818-01, 6.794095682990248-01,- 
1 2. 190863625159828-01, 4.333953941292478-01, 2.692667193099968-01,- 
1 1.488743389816318-01, 2.955242247147538-01/ 

DATA BB/- 

1 9.782286581460578-01, 5.566856711617408-02, 8.870625997680958-01,- 
1 1.255803694649058-01, 7.301520055740498-01, 1.862902109277348-01,- 
1 5.19096129206812B-01, 2 . 33 1 93764591990B-01 , 2.69543155952345E-01,- 
1 2.628045445102478-01,-0. 8-00, 2.729250867779018-01,- 

1 9.815606342467198-01, 4 .71753363865120B-02, 9.041172553704758-01,- 
1 1.069393259953188-01, 7.699026741943058-01, 1.600783285433468-01,- 
1 5.873179542866178-01, 2.031674267230668-01, 3. 67831498998 180E-0 1, - 
1 2. 334925365383558-01, 1.252334085114698-01, 2.491470458134038-01,- 
1 9.841830547185688-01, 4.048400476531608-02, 9.175983992229788-01,- 
1 9.212149983772808-02, 8.015780907333108-01, 1.388735102197878-01,- 
1 6.423493394403408-01, 1.781459807619468-01, 4.484927510364478-01,- 
1 2.078160475366898-01, 2.304563159551358-01, 2.262631802628978-01,- 
1-0. 8-00, 2.325515532300748-01, 9. 862838086968 12B-0 1 , - 

1 3.511946033175208-02, 9.284348836635748-01, 8.015808715976008-02,- 
1 8.272013150697658-01, 1.215185706879038-01, 6. 8729290481 16858-01 , - 
1 1.57203167158194E-01, 5.152486363581548-01, 1.855383974779388-01,- 
1 3.191123689278908-01, 2.051984637212968-01, 1.080549487073448-01,- 
1 2.152638534631588-01, 9.879925180204858-01, 3.075324199611708-02,- 
1 9. 372733924007C6E-01, 7.036604748810808-02/ 

DATA CC/- 

1 8.482065834104278-01, 1.071592204671728-01, 7.244177313601708-01,- 
1 1.395706779261548-01, S . 70972 1 726085398-0 1, 1.662692058169948-01,- 
1 3.941513470775638-01, 1.861610000155628-01, 2.011940939974358-01,- 
1 1.984314853271128-01,-0. 8-00, 2.025782419255618-01,- 

1 9.894009349916508-01, 2.715245941175408-02, 9.445750230732338-01,- 
1 6.225352393864808-02, 8.656312023678328-01, 9. 515851168249308-02,- 
1 7.554044083550038-01, 1.246289712555348-01, 6.178762444026448-01,- 
1 1.495959888165778-01, 4.580167776572278-01, 1. 69I56S19395003E-0 1 , - 
1 2.816035507792598-01, 1.826034150449248-01, 9.501250983763708-02,- 
1 1.894506104550688-01/ 

DO 50 I»1,70 
J«2*I 
JJ*2*I-1 
A«I) 

H(I)=C(J) 

50 CONTINUE 


QUA01 » 0. 


IP 

(11 

.EQ. 

XO) 

RETUBN 

N « 

UN 




IP 

( N 

.GT. 

16) 

N»16 

INDRTs 

BOO (N,2) ^ 

► 1 

INDKT 

= 1. 

H IS 

EVEN 

INDRT 

= 2. 

N IS 

000 

GO 

TO 

(20U, 

210) , 

,IHDNT 

BIN 

== (N 

•N)/« 

-1 


BAX 


•(»*2))/« 

-2 

60 

TO 

215 




210 BIN- (N*N-9)/4 ♦! 

WAX* (N* (N*2)-11)/4 
215 DELX-(XO-XI)/FLOAT(NS) 

HDELX*.5»DBLX 
DO 260 Rsl^tis 
XAaPLOAT (K-1)*DELX ♦ XI 
DO 220 1=* 1,2 
DO 220 J==BIN,HAX 
A(J)* -A(J) 

220 QOAD1 = QDAD1«B(J)*FCPX(XA>BDEI.X*(1.>A(J))> 


17 


IF(INDKT .Sy , 2) QDAD1 aQaADUH{MAX+1) *FOPX (X A*HOELX* ( 1. ♦ A (MAX + 1 ) ) ) 
260 CCKTINUE 

QUAD1 = HDELX*QUAD1 
270 RETURN 
END 
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APPENDIX D 


CALCULATION OF LIMITS OF INTEGRATION FOR NUMERICAL 
EVALUATION OF NORMALIZED SHAPE FUNCTION 

The integration indicated in equation (18) is performed numerically by replacing the 
infinite limits of integration and integrating the product Sg(i^, 0) over an interval 

where it is significant. Determination of practical limits involves finding an approxima- 
tion for i|/(x, 0) and for Sg(v, i/’)i//(x, 0). As these approximations involve products of 
Gaussian functions, it is necessary to find an expression for the product of two Gaussian 
functions. 

-a(x-x. -b(x -x, 

Consider the two Gaussian functions e ^ and e and their product 

exp |^-a(x - Xj)^J exp |^-b(x - = exp j^a(x - Xj^)^ + b(x - 

Expanding the exponent results in 

a(x - Xj^)^ + b(x - Xg)^ = a^x^ - 2xxj^ + x^^+ b^x^ - 2xxg + Xg^ 

= x^(a + b) - 2x(ax^ + bxg) + ax^ + bx| 


Completing the square in the first two terms gives 


(a + b) 


- ax, +bx« 

x2 - 2x _i 2 . + 

a + b 


(axj + 


bxg)^ 


(a + b)‘“ 


2,2 

+ axj + bxg - 


(axj + bxg)'’ 


a + b 


= (a + b) X - 


ax, + 


. ab 


+ (X, - x,)‘ 

a+b / a+b ^ ^ 


and 


exp j^-a(x - exp [^-b(x - X 2 )^j = exp 


ab 

a + b 


(Xi - Xg)^ 


exp 


-(a+b) X - 


axj^ + bXg) 


a + b 


(Dl) 
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The r^ht side of equation (Dl) is another Gaussian function, which has its maximum . 
value exp J- [ab/(a + bHCxj - at x = (axj + bx 2 )/(a + b) and has a bandwidth pro- 

portional to (a + b)"^/. 

The approximation of 0) involves replacing Sj^ of equation (8) with 


®R" 


J3 


n(v” - v' + 


exp 


4 In 2 /, tt 


L 4^2 


/ ” ’\2 
(v - v) 


= ejq) 


4 In 2 /.II 


/II ’\2 
(V - I/) 




(D2) 


The normalizing constants are discarded as only the functional behavior is now of 
concern. 

The Lorenz-type Raman distribution is replaced by a Gaussian distribution having 
the same bandwidth. Then the Voight function (eq. (9)) is aH>roximately 



Using the results of equation pi) gives 



This integral can be evaluated analytically. The result is 


'^^ .exp 


4 ln2 



P5) 


Equation p5) would indicate that the Voight bandwidth is given approximately by 



but a better approximation is given in reference 15. 


It is 




(D6) 
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Then 


Sy(p, ia) oc ejq) 


ii!il ipy-py 

4 


P7) 


and when equations (15) and (16) are used, 


i//(x, 0) exp 


ii2_i („’ _ p f 

2 •> 


L 


(D8) 


This approximation to the first Doppler broadening integral gives an approximation 
to the shape fimction. Substituting (D8) in equation (18) yields 



Again using equation (Dl) results in 



(DIO) 


At this point, it is possible to choose reasonable integration limits and use equa- 
tion (18) to evaluate S(v, ^j) numerically. Equation (DIO) indicates that the integrand of 
equation (18) behaves roi^hly like a Gaussian function centered at 


r2 *r2 
Ry + RS 


(Dll) 
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and with bandwidth 


+ — 

\4 4 

This information, together with the fact tiiat 


rl/2 


(D12) 





y^oo 


1 exp 

f,. or> 

_ 4 1n2 )2l 

L 4 J 

di» = 0. 999. . . i 

I exp 

r. 4 2 . „^,2l 

_ 4 


allows equation (18) to be replaced by 


dv 

(D13) 



v^+3Bw 


S(v, V,) = J- 


exp 


'C"3% 


4 In 2 

rI 


{v - vf 


0)di/ 


(D14) 


Equation (D14) is now in suitable form for integration by some standard technique, such 
as Simpson’s rule or Gaussian quadratures (ref. 16). The latter technique is chosen be- 
cause it gives the same results as a Simpson’s rule method that uses step size control, 
but is much faster. 
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Figure 1. - Line broadening. 




Figure 3. - lorentz broad^ing funrfior 

SR»(RR/2)/7r[(v- V())2 + (Rp/2)2]. 


24 
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Wavelength, m 

Figure 10. - Computed rotational Raman spectrum for N 2 at 1100 K and 1 atmosphere. 



Wavelength, m 


Figure 11. - Computed and measured Raman spectra for No at 300 K and 19.4 atmospheres, ftr both curves, laser 
bandwidth. 0.166 centimeter’*; siit bandwidth, 3. 31 centimeter'*; recipraal of wave number of laser line, 
5145xl0’*0 meter (5145. A). 
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